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ABSTRACT 



Context. The evolution of radial gradients of metallicity in disk galaxies and its relation with the disk formation are not well un- 
derstood. Theoretical models of galactic chemical evolution make contrasting predictions about the time evolution of metallicity 
gradients. 

Aims. To test chemical evolution models and trace the star formation and accretion history of low luminosity disk galaxies we focus 
on the Local Group galaxy M33. 

Methods. We analyze O/H and S/H abundances in planetary nebulae, Hn regions, and young stars, together with known [Fe/H] abun- 
dances in the old stellar population of M33. With a theoretical model, we follow the time evolution of gas (diffuse and condensed in 
clouds), stars, and chemical abundances in the disk of M33, assuming that the galaxy is accreting gas from an external reservoir. 
Results. Our model is able to reproduce the available observational constraints on the distribution of gas and stars in M33 and to 
predict the time evolution of several chemical abundances. In particular, we find that a model characterized by a continuous infall of 
gas on the disk, at a rate of M m f w 1 M Q yr , almost constant with time, can also account for the relatively high rate of star formation 
and for the shallow chemical gradients. 

Conclusions. Supported by a large sample of high resolution observations for this nearby galaxy, we conclude that the metallicity in 
the disk of M33 has increased with time at all radii, with a continuous flattening of the gradient over the last ~ 8 Gyr. 

Key words. Galaxies: abundances, evolution - Galaxies, individual: M33 



1. Introduction 

Local Group (LG) galaxies, with their variety of morphological 
types, are ideal candidates for testing the predictions of galactic 
chemical evolution (GCE) models. The advent of a new gener- 
ation of telescopes and sensitive instruments, such as wide field 
cameras and high resolution spectrographs, has given detailed 
information about the distribution of atomic and molecular gas, 
stellar content, chemical abundances from Hn regions, planetary 
nebulae (PNe), and stars in external galaxies, especially in the 
LG. 

M33 (NGC 598) is a low-luminosity late-type spiral galaxy 
(morphological type Sc 11-111, cf. van den Bergh 2000) located 
in the LG at a distance of 840 kpc (Freedman et al. |1991| l. Its 
proximity implies a large angular size which, together with its 
modest inclination (i = 54°), makes this galaxy particularly suit- 
able for detailed studies of its stellar content, gas distribution, 
and abundance gradients. 

M33 is an ideal target to test GCE models and constrain 
galaxy assembly processes because it shows no signs of recent 
mergers and no presence of a prominent bulge or bar compo- 
nent (Regan & Vo gel [19941 1. If mergers were a recent dominant 
process then an extended stellar halo with several substructures 
should be visible (Bullock & Johnston 2005). Spectroscopic sur- 
veys of red giant branch (RGB) stars in M33 show no clear ev- 
idence for these substructures (Ferguson et al. 2006), indicat- 
ing that M33 is evolving relatively undisturbed. In addition, star 
counts in the outer disk and individual stellar velocity disper- 
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sion measures in the optical disk indicate that any stellar halo 
component contributes for less than a few percent to the total 
luminosity (Barker et al. 2006, McConnachie et al. 120061) . 

The recent identification of individual stars in M33 (e.g. 
Block et al. 120061 Barker et al. 120061 ) suggests the existence of 
different stellar populations and therefore of separate episodes 
of star formation, but a satisfactory understanding of the overall 
evolution of M33 is still lacking. Chemical abundances of stellar 
populations of different ages, fundamental constraints for GCE 
models, are becoming available for this galaxy (e.g. Magrini 
et al. 12TJ041 120071 Urbaneja et al. 120051 Crockett et al. 120061 
Barker et al. 2006). Surveys of molecular hydrogen through the 
CO J=l-0 line (Engargiola et al. |2003l Heyer et al. 120041 ), to- 
gether with limits on stellar masses from dynamical mass model- 
ing (Corbelli 2003 ) and photometric surveys, now allow a quan- 
titative comparison between observable quantities and model 
predictions. 

Early attempts to model the chemical evolution of M33 date 
back to the seminal work of Diaz & Tosi ( 1984 ), and were later 
developed by Molla et al. ( 1996 ) through a formalism originally 
devel oped b y Ferrini et al. d!992| l. Molla et al. d 19961 ) and Molla 
et al. ( 1 19971 ) tackled the issue of the time evolution of radial dis- 
tributions of elemental abundances in disk galaxies of different 
morphological type, including M33, introducing in their model 
a radial dependence of the infall rate and the star formation rate 
(hereafter SFR). A similar approach was also taken by Goetz 
& Koeppen ( 119921 ), and Koeppen (1199411 wi th analytical GCE 
models. More recently, Molla & Diaz (2005) parameterized in- 
puts and observational data relative to disk galaxies in order to 
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compute a grid of multiphase chemical evolution models able 
to predict the evolution of irregular and spiral galaxies. Their 
predictions compare well (in a statistical sense) with the avail- 
able data for several nearby galaxies, including M33. The mod- 
els for M33 quoted above need however to be reconsidered at the 
light of the new data that have recently become available, such 
as radial gas profiles and the metallicity determinations in stellar 
populations of different ages (e.g. Barker et al. 2006, Magrini et 
al. 120041 1. The goal of this paper is therefore to model in detail 
the evolution of M33 in order to reach a satisfactory agreement 
with the full set of available observations. In particular, we wish 
to investigate the controversial issue of the time evolution of ra- 
dial metallicity gradients. 

The problem of the time variation of radial chemical gra- 
dients is far from being settled, either theoretically or observa- 
tionally, even in the case of our own Galaxy (see e.g. Goetz & 
Koeppen [T9921 Ko eppen [19941 T osi [19961 Molla et al. [19971 
Henry & Worfhe y [T^99l Maciel 120001 Maciel et al. [20031 
Stanghellini et al. 2006). Different models predict opposite be- 
haviours of the metallicity gradient, showing the sensitivity of 
this observable to the adopted parameterization of the physical 
processes. In this respect, GCE models can be divided in two 
groups: those where the metallicity gradients steepen with time 
(Tosi [19881 Chiappini et al. [19971 Chiappin i et al [200TT) and 
those where they flatten with time (Molla et al. 119971 Portinari & 
Chiosi [19991 Hou et al.[2000). The main difference between the 
two groups is how fast chemical enrichment proceeds in inner 
and outer regions of the galactic disk. In models of the former 
group, the outer disk is pre-enriched by the halo, and its metal- 
licity is affected very little by the subsequent low star formation 
activity. In contrast, in models of the latter group, a vigorous star 
formation activity in the inner disk results in a rapid increase 
of the metallicity near the galactic center. Here the metallicity 
increases very fast reaching its final value in the first 2-3 Gyr 
of disk evolution, whereas the outer disk enrichment increases 
much slowly and progressively flattens the radial gradient. 

Observationally, it is possible to discriminate between these 
two groups of models by comparing the metallicity gradient 
of the "old" stellar population (RGB stars and PNe) with the 
present-day gradient outlined by Hn regions and young stars. In 
particular, the comparison of PNe and Hn regions investigates 
the similarities between the observational techniques and abun- 
dance analysis for the two samples (Maciel et al. 2005, 2006). 
In our Galaxy the time variation of the gradients is mainly de- 
rived from chemical abundances of PNe, open clusters, Cepheids 
and young objects, such as OB associations of stars and Hn re- 
gi ons (e. g. Maciel et al. [20031 [20031 [20061 Friel [19951 Friel et 
al. 120021 ). Observationally, the [Fe/H] gradient in open clusters 
supports a time flattening of the abundance gradient in agree- 
ment with the gradients derived from Cepheids of well deter- 
mined age. According to Maciel et al. (2005 ), the average flatten- 
ing rate for the last 8 Gyr is about 0.005-0.010 dex kpc'Gyr 1 . 

In M33, it is not possible to derive the time evolution of the 
gradient from abundances of PNe alone, as in the case of the our 
own Galaxy, because of the difficulty to sub-divide PNe in differ- 
ent age groups (Magrini et al. 2004). However, a comparison be- 
tween abundances measured in Hn regions (tracing present-day 
abundances in the interstellar medium, ISM) and in PNe (tracing 
the ISM composition at the time of formation of their progeni- 
tors) should in principle reveal any change in the radial distri- 
bution of the elements occurred in the last few Gyr, for those 
elements not affected by low-mass stellar nucleosynthesis like 
O, S and Ne. 



The paper is organized as follows: in Sect. [2] we collect and 
critically examine the observational data adopted to constrain 
our GCE model described in Sect. [3j our results are presented 
in Sect. [4] with a particular attention to the issue of the time 
evolution of the metallicity gradients (Sect. [5]); finally, in Sect. [6] 
we summarize our conclusions. 



2. Observational data 

In this section we examine the observational constraints for the 
evolutionary model of M33, namely: the radial profiles of gas 
and stars, the SFR, and the chemical abundances. 



2. 1 . Gas surface density 

The gas and stellar mass surface densities in M33 have been an- 
alyzed by Corbelli (2003) using the gas kinematics. The neu- 
tral hydrogen surface density profile has been derived from high 
sensitivity observations with the Arecibo-305m telescope after 
applying corrections for the disk inclination with respect to the 
line-of-sight according to the tilted ring model described in de- 
tail by Corbelli & Salucci d2000l) . 

Imaging of molecular clouds complexes has recently been 
carried out in M33 with the BIMA interferometer and with 
the FCRAO-14m telescope using the CO J=l-0 line transition 
(Corbelli l2003l Engargiola et al. 120031 Heyer et al. 120041 . These 
data sets reveal the presence of both compact, massive molecular 
clouds, and of a more distributed molecular gas component. The 
total molecular mass is ~ 2 x 10 8 M Q , corresponding to ~ 10% 
of the ISM mass of this galaxy. Giant molecular clouds (GMCs) 
detected by the interferometer account for only ~ 20% of the 
molecular gas mass. Most likely some of the total molecular 
mass detected by FCRAO just resides around the complexes im- 
aged by the interferometer while some is in diffuse form spread 
out throughout a more extended area. 

The azimuthally averaged surface density of atomic and 
molecular gas, corrected for the disk inclination with respect to 
the line-of-sight, from data in Corbelli (2003), can be well rep- 
resented as follows: 



-gas 



21 exp 



R 



7.8 kpc 



M pc~ 



(1) 



2.2. Stellar surface density 

The main stellar component of M33 is a thin disk: the central 
bulge and the stellar halo contain less than 10% of the total lu- 
minous mass. A'-band photometry, which more closely reflects 
the underlying stellar mass distribution in the disk being less af- 
fected by extinction, points out to an exponential stellar light dis- 
tribution declining radially with a scale length of 1 .4 kpc (Regan 
& Vogel ll994l ). Both B- and /T-band photometry show that de- 
partures of the stellar light profile from an exponential law are 
visible in the innermost kpc, where either the disk scale length 
changes, or a compact bulge is present. Using the Zf-band scale 
length, together with data on the atomic and molecular gas dis- 
tribution, Corbelli (2003 ) has inferred the stellar mass-to-light 
ratio, and therefore the stellar mass density of this galaxy. This 
includes also the mass in stellar remnants, being derived via a 
dynamical mass model. The best fitting three component model 
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(gas, stars, dark matter) to the rotation curve requires a stellar 
surface mass density of 



Sstars = 430 exp 



R 



1.42 kpc 



M pc 



-2 



(2) 



The rotation curve of M33 rises in the center out to a galacto- 
centric radius R w 3 kpc, then it flattens out up to the outermost 
measured point. Radial mixing flows, possibly driven by shear 
generated by differential rotation, can be present at R > 3 kpc. 

2.3. The SFR 

The blue colour of M33 (B — V = 0.55) and the prominent Hn 
regions indicate vigorous star formation activity in the nuclear 
region and along the spiral arms, where the largest molecular 
complexes are found (Engargiola et al. 2003). In general, optical 
and infrared luminosities can be used to trace the SFR in external 
galaxies. Optical observations rely on the conversion of Ha lu- 
minosities into SFRs (Kennicutt 1998). The big uncertainty with 
this method is the Ha luminosity correction for dust extinction. 
Other methods involving the far infrared (FIR) luminosity are 
less affected by extinction. In this case, however, the conversion 
factor for deriving the SFR depends sensitively on the assump- 
tion about the fraction of the total luminosity re-irradiated in the 
infrared (IR) and also on the type of stars contributing to dust 
heating. The end result is that the IR conversion factor is sub- 
ject to large uncertainties. Fortunately, M33 is known to have 
moderate extinction and the optical method seems appropriate 
for this galaxy. Extinction for M33 has been studied by Israel & 
Kennicutt ( 1980) using a sample of bright Hn regions and free- 
free radio emission measurements. They found a clear radial de- 
pendence of the extinction law. Devereux et al. ([1997) on a larger 
sample of Hn region, using the same radio-Ha method, found no 
clear radial trend but an average extinction of A v m 1 mag. We 
will use this value without introducing any other "uncertain" ra- 
dial dependence. 

We determine the Ha surface brightness profile from the Ha 
image of Hoopes & Walterbos (2000), averaging the Ha flux 
along ellipses. These are projections of circular rings 0.5 kpc 
wide, inclined by 54° at a position angle 23°. The diffuse ion- 
ized gas (DIG) contributes for ~ 40% to the total observed Ha 
flux (Hoopes & Walterbos 2000) and we assume that the DIG 
emission is unaffected by dust extinction. DIG emission is how- 
ever important in computing the SFR since field O and B stars 
are responsible for it. After subtracting 5% of the recovered flux 
due to possible [Nn] lines contamination, and correcting 60% 
of the emission for 0.83 mag of extinction (a factor 1.2 smaller 
than extinction in the visual), the surface brightness is converted 
in SFR density following Kennicutt et al. ( 1998 ), 



SFR 



L(Ha) 



1.26 x 10 41 ergs" 



M Q yr 



(3) 



The resulting SFR per unit area is shown in Fig. [2] The SFR 
integrated over the whole disk is ~ 0.4 M Q yr , comparable to 
the value given by Kennicutt ( 1998 ) and also to the value derived 
from the FIR luminosities (see Hippelein et al. 2003 and refer- 
ences therein). We shall use the SFR derived by Kennicutt ( 1989 ) 
from Ha luminosities of bright Hn complexes and the FIR SFR 
by Heyer at al.(2004) respectively as lower and upper limits for 
our model results. 

The SFR in the range from 1 to 8 Gyr ago can be es- 
timated from the PNe population as described by Magrini et 
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Fig. 1. The time evolution of the infall rate for the accretion 
model (panel a) and for the collapse model (panel b). Curves 
represents the infall rate at 0.5 {dotted curve), 2 (dashed curve), 
3 (long-dashed curve), 5 (dot-dashed curve), 8 (long dash-dotted 
curve), 12 (long-short dashed curve), and at 13.6 Gyr (solid 
curve). 



al. d2005al l. In short, the method consists in counting PNe, es- 
timating the mean mass of the PNe progenitors and obtaining a 
measure of the total mass of the intermediate-age stellar popula- 
tion. As a result, the intermediate-age stellar mass is proportional 
to the number of PNe. Extrapolating the number of known PNe 
(Ciardullo et al. 120041) in M33 within 4 m ag from the cut-off 
of their luminosity function (Ciardullo et al. 1 1989ft , and consid- 
ering a time interval of 3 Gyr for the formation of their pro- 
genitors (the youngest, since we are considering the brightest 
PNe) we derive a mean SFR=0.55 M Q yr _1 ,(l^l Gyr ago) inte- 
grated over the whole disk. On the other hand, considering older 
PNe, extrapolating within 8 mag from the cut-off of their stan- 
dard luminosity function, and considering the whole time inter- 
val for the formation of progenitor stars, we obtain a higher mean 
SFR=1.1 M Q yr 1 (1-8 Gyr ago). 

2.4. The infall rate 

Star formation in M33 extends much beyond the region where 
molecular complexes are located. Carbon stars have been found 
even at larger radii where the drop of the Ha flux occurs (R « 
6 kpc, Kennicutt |1989l l. The presence of such an intermediate- 
age stellar population in the outer disk suggests that accretion 
of gas is still taking place (Davidge 2003; Block et al. 2004 
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Fig. 2. Comparison of the predicted time evolution star forma- 
tion rate (line types as in FigJTJ with observations for the ac- 
cretion model (panel a) and with the collapse model (panel 
b). Data are from Kennicutt et al. d 19891 crosses); Hoopes & 
Walterbos (120001 filled circles); Heyer et al. (120041 triangles). 



2006; Rowe et al. 2005 ) and that the disk of the galaxy is form- 
ing inside-out. A likely source of the accreted gas could be the 
cosmic web (Keres et al. 2005 ), the warm/hot filaments extend- 
ing between galaxies, where most of the baryons reside. 

Is there any evidence that M33 still in the process of accret- 
ing mass from the intergalactic medium? A conspicuous cloud of 
neutral gas infalling into the disk of M33 recently detected at 21- 
cm by Westmeier et al. (2005) supports this picture. The cloud 
might have been tidally stripped from a distant faint dwarf or 
dark galaxy (with no associated stellar counterpart) or it might 
just be gas condensing out of intergalactic medium filaments. 
This finding confirms that M33 is still in the process of accret- 
ing gas to fuel star formation. The gas infall rate over a disk re- 
gion corresponding to the cloud size is about 2 M pc~ 2 Gyr -1 . 
There are however observational difficulties in estimating the to- 
tal amount of gas infalling into the M33 disk. Velocities of the 
infalling clouds are expected to be similar to velocities of high 
velocity clouds related to the MW disk, and therefore confusion 
limits our ability to distinguish the two different populations in 
the small sky area where M33 is located. Any high velocity cloud 
in projected proximity to the M33 disk, placed at the distance of 
M33 would result in a much higher infall rate over M33 than 
what could be inferred for the same cloud if it were closer to us 
and infalling into the Milky Way (MW) disk. However, the num- 
ber of high velocity clouds throughout the whole sky is large, 



and most of them are connected with gas infalling into the MW 
galaxy. Observational estimates of the infall rates into the MW 
disk range between 0.1 and 5 M G yr~' (e.g. Wakker et al. 1999; 
de Boer 2004). Rates « 1 M yr~' or larger are favoured by con- 
siderations on the rapid molecular gas depletion, on the G-dwarf 
problem, and by some recent evolution models of disc galaxies 
(e.g. Naab & Ostriker 2006 and references therein). The total in- 
fall rate for M33 is expected to be lower than in the MW because 
its total mass is lower than the mass of the MW. This consider- 
ation applies both to halo gas infall models, which involve the 
collapse time scale, and to models where the gas starts flowing 
in from the intergalactic medium (see Section 3.1). Given that 
quantitative observations of the infall rate over a wide area of 
M33 and a wide interval of time are difficult to obtain, in Sect. [4] 
we shall discuss the implications of models characterized by dif- 
ferent infall rates. 



2.5. Chemical abundances: stars and ISM 

Aller (1 19421 ) first obtained spectra of Hn regions in M33 and 
derived a radial gradient of the [On]/[Om] line ratio, that he at- 
tributed to an excitation gradient. Smith (1975 ), in subsequent 
spectroscopic study, found clear evidence for a metallicity gra- 
dient. Further spectroscopic studies of Hn regions were carried 
on by Kwitt er & Aller (fTMTT l and Vilchez et al. (fT988b . Garnett 
et al. (119971 ) using data from previous observations with derived 
electron temperatures, re-determined the abundances obtaining 
an overall O/H gradient of -0.1 1 + 0.02 dex kpc 1 , including the 
central regions of M33. The steepness of the gradient is how- 
ever dominated by the innermost data point. Subsequent stud- 
ies of the metallicity gradient which do not include the central 
regions obtain a much shallower gradient. Recently, Willner & 
Nelson-Patel (2002) have derived Ne/H abundances for 25 Hn 
regions from infrared lines, obtaining a [Ne/H] gradient -0.034+ 
0.015 dex kpc -1 , outside -0.5 kpc from the center. Crockett et 
al. ([2006), measuring chemical abundances in six Hn regions, 
have derived a [Ne/H] gradient of -0.016 + 0.017 dex kpc 1 and 
a [O/H] gradient of -0.012 + 0.011 dex kpc 1 , much shallower 
than in previous studies, but at the limit of significance due to the 
smallness of the sample. All these determinations are consistent 
with each other and imply that excluding the innermost ~0.5 kpc 
region where a bulge might be present, the metallicity gradient 
of Hn regions in the disk of M33 is very shallow. 

Chemical abundances of PNe have been derived by Magrini 
et al. ( 120041 ) via optical spectroscopy and photoionization mod- 
eling. Stellar abundances have been obtained by Herrero et 
al. d 19941 ) for A B-type supergiant stars, McCarthy et al. (119951 ) 
and Venn et al. ( 1 19981 ) for A-type supergiant stars, Monteverde 
et al. (fT997l 12000] ) and Urbaneja et al. ( 120051 ) f or B-ty pe super- 
giant stars. The larger sample by Urbaneja et al. (2005 ) indicates 
an [O/H] gradient of -0.06 + 0.02 dex kpc 1 . Very recently the 
detection of beat Cepheids, which are young intermediate mass 
stars (Beaulieu et al. 2006 ) allowed the derivation of their metal- 
licity making use of stellar pulsation models and of the mass- 
luminosity relation. Beaulieu et al. (2006) derived a metallic- 
ity gradient of -0.20 dex kpc 1 , and inferred an [O/H] gradient 
-0.16 dex kpc 1 . Metallicities of the old stellar population have 
been obtained via deep CCD photometry and colour-magnitude 
diagrams by Stephens & Frogel (20021, Kim et al. (120021 ) in 
inner fields, and by Galletti et al. (|2"004jl, Tiede et al. d2"004j l. 
Brooks et al. (2004) and Barker et al. (|2006 ) in outer fields. The 
[Fe/H] gradient related to the whole RGB stellar population has 
a slope -0.07 + 0.01 dex kpc 1 (Barker et al. |2006l 



Magrini et al.: The evolution of M33 



5 



3. The model 

The GCE model adopted in this work is a generalization of 
the multi-phase model by Ferrini et al. (11992ft , built for the so- 
lar neighborhood, and subsequently extended to the entire MW 
(Ferrini et al. [1994), and to other disk galaxies (e.g. Molla et 
al. [19961 [19971 Molla & Diaz 120051 Here we summarize the 
main characteristics of the model (see Ferrini et al. [1992 for de- 
tails). 

The galaxy is divided into N coaxial cylindrical annuli with 
inner and outer galactocentric radii (i = l,N) and 
respectively, mean radius /?,+i/2 = (R; + Ri+\)/2 and height 
h(Ri + i/2). Each annulus is divided into two zones, the halo and 
the disk, made of diffuse gas g, clouds c, stars s and stellar rem- 
nants r. The halo component is intended here quite generally as 
the primordial baryonic halo, or as material accreted from inter- 
actions with small LG galaxies or from the intergalactic medium 
during the life-time of the galaxy. In the following, a subscript H 
or D indicates the halo or the disk, respectively. For a spherical 
halo of radius Rn+\ = Rh, 



h(R i+m ) = (R 2 H -Rl U2 ) l/2 , 

and the halo volume in the i— th annulus is then 

V H j = n(R 2 M -R 2 )h(R M i2). 



(4) 



(5) 



If zd is the disk scale height (assumed independent on radius), 
the volume of the disk in each annulus is 



V DJ = n{R 2 M - Rf)z D . 



(6) 



At time t = all the mass of the galaxy is in the form of diffuse 
gas in the halo. At later times, the mass fraction in the various 
components is modified by several conversion processes: diffuse 
gas is converted into clouds, clouds collapse to form stars and are 
disrupted by massive stars, stars evolve into remnants and return 
a fraction of their mass to the diffuse gas. The disk of mass Mo{t) 
is formed by continuous infall from the halo of mass M//(?) at a 
rate 



dM L 
dt 



: fM H , 



(7) 



where / is a coefficient of the order of the inverse of the infall 
time scale. Clouds condense out of diffuse gas at a rate p. and are 
disrupted by cloud-cloud collisions at a rate H' , 



dM c 
~dT 



= fiM 3 ' 2 - H'M 2 , 



(8) 



where M g (t) and M c (t) are the mass fractions of diffuse gas and 
clouds, respectively. Stars form in the halo and the disk by cloud- 
cloud collisions at a rate H and by the interactions of massive 
stars with clouds at a rate a, 



dM s , 

= HM 2 C + aM s Mc - DM S 

dt 



(9) 



where M s (t) is the mass fraction in stars and D is the stellar death 
rate. 

Each annulus is then evolved independently (i.e. without ra- 
dial mass flows) keeping fixed its total mass from t = to 
fgai = 13.6 Gyr computing the fraction of mass in each com- 
ponent in the two zones, and the chemical composition of the 
gas (assumed identical for diffuse gas and clouds). The rate co- 
efficients of the model are all assumed to be independent on time 
but functions of the galactocentric radius R. Their radial depen- 
dence for a model of the MW is discussed in detail by Ferrini et 



Distance 


s4U kpc 


Type 


C~ TT TTT 
oC 11-111 


Tnpl in i nn 

lll^lllltlllVMl 


54° 


Position angle 


23° 


Luminosity 


6.5 x 10 9 L G 


Maximum rotation speed 


130 km s- 1 


Optical radius R opt 


6.6 kpc 


Equivalent Solar radius 


3.9 kpc 



Table 1. M33: adopted physical characteristics 



al. ( fl994l l, Molla et al. ( 1996), and Molla & Diaz (2005). In gen- 
eral, coefficients representing condensation processes (like e.g. 
the formation of clouds from diffuse gas), being proportional to 
the inverse of the dynamical time, scale with the inverse square 
root of the zone volume, whereas the coefficients of binary pro- 
cesses (like e.g. the formation of stars by cloud-cloud collisions) 
scale as the inverse of the zone volume; the coefficient of star 
formation induced by stars is independent on radius. 

The model described in this section can be applied in gen- 
eral to study the chemical evolution of any disk galaxy. We now 
discuss the values of the various coefficients that we have specif- 
ically adopted for M33. As a starting point, we have adopted the 
set of values by Ferrini et al. (119921 ) for their "best model" of 
the solar neighborhood, applying them at the equivalent solar ra- 
dius Q The values at any other radius are then determined by the 
scaling relations described in Sect. [3] 

The mass fractions in each component computed by the 
model in each annulus are then converted into surface densities 
using as a normalization the total surface density, (i.e. the sum 
of the Hi, H2, stellar and remnants surface densities, as given by 
equations Q] and [2]). In the comparison with the data, we iden- 
tify the gas and cloud components with the Hi and the H2 gas, 
respectively. 

3.1. The infall coefficient 

We assume an exponential radial dependence of the infall rate /, 
f(R) = f exp(-R/A D ) (10) 

where fo is the infall rate at R = and Ad is a typical disk 
scale length (cf. Ferrini et al. 119941 Portinari & Chiosi I19991 ). 
In principle, Ad can be set equal to the scale length of the disk 
surface brightness in the B-band (~ 1.9 kpc, Freeman 119701 ), 
or in the K-bund (~ 1.4 kpc, Regan & Vogel[[994), or to the 
CO scale length (~ 2.7 kpc, Corbelli 120031 . In our model, the 
value of Ad that better reproduces the observations (gas and stars 
surface density profiles, radial SFR, see Figs. [2] [3] [4] [5]) is Ad ~ 
2.0 kpc. Larger values of Ad imply a huge amount of neutral 
hydrogen in the central regions, while lower values depress the 
SFR in the inner regions. 

The infall coefficient fo is usually set of the order of the in- 
verse of the collapse time t c for the galaxy, where t c depends 
on the total mass of the galaxy Mc a i as t c oc M g ^ 2 (Gallagher 
et al. 119841 ). This choice corresponds to a scenario where the 
disk is formed very rapidly at the beginning of a galaxy's life 
and we shall refer to this model as the collapse model. Since 
the mass of M33 is about 0.2 times the mass of the MW, the 
value of /o for M33 should be a factor ~ 2.2 smaller than in 
the case of the MW. Scaling the infall value adopted by Ferrini 
et al. ( [19921 ) for the MW (/ = 0.7 see Table [2]), corresponding 
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to a dynamical collapse time of the order of 10 yr at the solar 
radius, we obtain for M33 /o ~ 0.3. Over a disk of 10 kpc in ra- 
dius this choice results in a present-day average accretion rate of 
Mi„f « 0.05 M Q pc~ 2 Gyr 1 (FigQJ)). Thus our collapse model 
is intended to reproduce the characteristics of a rapid phase of 
disk formation, following the dynamical collapse of an extended 
halo. 

However, estimates of /o for the MW based on chemical con- 
straints (G-dwarfs metallicity distribution, O/Fe vs. Fe/H rela- 
tions, etc.) suggest lower values, /o « 0.012 for the MW, corre- 
sponding to a time scale of the order of several Gyr (see Molla 
& Diaz 2005 and references therein). Lower values of /o implies 
higher infall rates at the present time. If this is the case for M33 
it is unlikely however that the halo is the reservoir of the gas in- 
falling into the disk because there are no observational evidences 
for a gaseous halo. A possibility for supporting a high infall rate 
at present time, and in general an infall rate almost constant in 
time, is to consider a slow continuous accretion process from 
the environment. This is obtained for example using fit as 0.003, 
corresponding to an infall rate of 1.2 M yr~' in the M33 disk at 
present time or about and average of M^f « 3.8 M pc~ 2 Gyr -1 
(Fig [It)- We shall refer to this model as the accretion model. 

High infall rates are favoured by some recent papers describ- 
ing the evolution of disk galaxies (e.g. Naab & Ostriker 2006). 
Infall rate in the MW for example are predicted to be of order 2— 
4 M yr , with little variations over the last 10 Gyr. These val- 
ues are in good agreement with those inferred from observations 
of high velocity clouds and explain specific abundance patterns 
such as the high deuterium abundance at the solar neighborhood 
and at the Galactic Center (Chiappini et al. 2001). In these mod- 
els the intergalactic medium is often identified as the source of 
the gas infalling into the disk. Hydrodynamical simulations of 
galaxy formation and evolution predict that accretion of 'cold' 
gas from intergalactic filaments is present and becomes the dom- 
inant process for low mass galaxies and in low density regions. 
Galaxies of mass similar to M33 should accrete gas at a rate 
which decreases slowly with time and is currently » 1 M Q yr _I 
(Keres et al. 2005). This finding supports the infall rates used 
in our accretion model. Higher mass halos have larger accre- 
tion rates but a large fraction of the accreted gas has been shock 
heated at the time the galaxy formed, and does not flow in from 
large distances, as in the case of intergalactic filaments feeding 
the disk. The intergalactic feeding mechanism is well suited for 
M33 because of the absence of a massive gaseous and stellar 
halo, relic of the galaxy formation epoch and driver of disk ac- 
cretion processes. 

3.2. The star and cloud formation efficiencies 

The parameters regulating the conversion of diffuse gas into 
clouds and clouds into stars are: p (cloud formation from diffuse 
gas), H (star formation from cloud collisions), and H' (cloud 
dispersion). The efficiencies of these processes adopted in our 
models, together with previous values adopted for the MW and 
M33, are shown in Table [2] These efficiencies are in general 
functions of the morphological type of the galaxy (cf. Sandage 
[19861 Gallagher et al. [I9g4l : Ferrini & Galli JlW8t analyzed 
the behaviour of these parameters in the different morphological 
types of spiral galaxies, and found a reduction of p and H from 
a Sbc galaxy, like the MW, to a Scd galaxy, like M33. These 
scaling relations can be also compared with those of Molla & 
Diaz (2005) who studied spiral galaxies of different morpho- 
logical types. For a spiral galaxy of morphological type 6, such 
as M33, they adopted the following values (see their Table 2): 



en - 0.01 and = 0.15, which translate, at the equivalent solar 
radius, into H = 0.012 and p = 0.017 in units of 10~ 7 yr , us- 
ing the relations between e#, and H and between £ M and p given 
by Ferrini et al. (|1994| l. The H parameter adopted in the present 
work, in particular for the accretion model, has a higher value 
than t hat used for M33 by Molla & Diaz (|2003T > and by Molla 
et al. dl996|199Tb . This is probably due to the different param- 
eterization of the IMF used: Molla and collaborators used the 
Ferrini 's et al. (1 19901 1 IMF which predicts a larger number of 
massive stars with respect to other IMF, such as the Kroupa's et 
al. ([T993 ) IMF. Therefore their model needed a lower H value in 
order to reproduce the observed metal content. 

Note that in order to reproduce the observed lower content 
of molecular gas (Fig. [5} of M33 with respect to that found in 
Sc galaxies, a very low cloud formation efficiency is needed. 
The low molecular fraction in this galaxy is in agreement with 
the tendency of the gas fraction in molecular form to decrease 
from Sc to Irr morphological types, despite the higher SFR of 
the latter. 

3.3. The cloud dispersal coefficient 

The cloud dispersal coefficient H' is a measure of the probability 
that collisions between gas clouds result in the disruption of the 
clouds themselves and return cloud material to the diffuse phase. 
In principle, the rate of cloud dispersal should depend on radius 
in the same way as the rate of star formation by cloud-cloud 
collisions, namely as the inverse of the disk volume Vq (see 
Ferrini et al. 119941 ). However, this dependence makes impossible 
to reproduce the observed distributions of Hi and H2, because at 
small galactocentric radii, where the factor 1/Vd enhances the 
rate of formation of diffuse gas by cloud collisions, the Hi gas is 
underabundant (see Fig. |4j. Our model requires instead that the 
coefficient of cloud dispersal, H', should be roughly independent 
on radius. Since the dependence of the frequency of collisions 
from the inverse of the volume is a general characteristic of any 
binary collision process, an additional hypothesis must be made 
about the efficiency of the process of cloud dispersal. Namely, 
our modeling requires that the clouds in the outer galaxy are 
more efficiently dispersed by collisions than clouds in the inner 
galaxy. 

This model requirement is related to the observed absence 
of GMCs at large radii where star formation is still taking place, 
supporting the claim that clouds are easily destroyed in the outer 
disk of M33. This may be due to a radial change of some in- 
trinsic property of the clouds (e.g. compactness) or to additional 
processes that regulate the molecular/atomic hydrogen ratio in 
galaxies of low molecular fraction such as M33, like photodis- 
sociation by a pervasive interstellar radiation field (Elmegreen 
119931 Heyer et al.|2004). Because of the larger photon mean free 
path at larger radii, this process probably dominates the conver- 
sion of molecular clouds into atomic diffuse gas in the outer disk 
ofM33. 

3.4. Stellar yields 

We model the chemical enrichment of the gas using the formal- 
ism developed by Talbot & Arnett dl9731 l who introduced the 
restitution matrices Q/j(M,Z). The elements of these matrices 
are defined as the fraction of the mass of an element j initially 
present in a star of mass M and metallicity Z that it is converted 
in an element i and ejected. At each stellar mass M and metallic- 
ity Z corresponds one matrix Qij(M,Z). Our GCE model takes 
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into account two different metallicities, Z = 0.02 and Z = 0.006, 
and 22 stellar masses (21 for Z = 0.006), ranging from 0.8 to 
100 M©, for a total of 43 restitution matrices. 

We update the nucleosynthesis yields used to build the 
Qij(M,Z) matrices in the following way. For low- and 
intermediate-mass stars (M < 8 M ) we use the yields by 
Gavilan et al. (2005 ) for both values of the metallicity. We have 
also tried the yields of Marigo (2001 ), which however do not 
give any appreciable difference in the computed gradients of 
chemical elements produced by intermediate mass stars, as N. 
For stars in the mass range 8 M < M < 35 M we adopt the 
yields by Chieffi & Limongi (2004) for Z = 0.006 and Z = 0.02. 
The yields of more massive stars are affected by considerable un- 
certainties associated to different assumptions about the model- 
ing of processes like convection, semiconvection, overshooting, 
mass loss. Other difficulties arise from the simulation of the su- 
pernova explosion and the possible fallback after the explosion, 
that strongly influences the production of iron-peak elements. It 
is not surprising then that the results of different authors (e.g. 
Arnett [T9931 Woos ley & Weaver [l995l Thielemann et al. [T996l 
Aubert et al. 1996) disagree in some cases by orders of magni- 
tudeQ. In our models, we estimate the yields of stars in the mass 
range 35 M© < M < 100 M by linear extrapolation of the yields 
in the mass range 8 M Q < M < 35 M Q . 

3.5. The IMF 

Recent work supports the idea that the IMF is universal 
in space and constant in time (Wyse 1997} Scalo [1998; 
Kroupa 2002), apart from local fluctuations. As reviewed by 
Romano et al. (120051 there are several parameterizations of 
the IMF that have been used in GCE models, starting from 
the original Salpeter ( 1 19551 1 power-law: piecewise power-laws 
(Tinsley [19801 Scalo [19861 Kroupa et al. [19931 Scalo [19981 . 
polynomial approximations (Ferrini et al. 1990), and logarithmic 
polynomials (Chabrier 2003). Romano et al. ( 120051 ) studied the 
sensitivity of GCE models to different parameterizations of the 
IMF and found that the Scalo ( [19861 ), Kroupa et al. ([19931 and 
Chabrier (2003 ) IMFs are generally more consistent with obser- 
vational data than those of Salpeter ( 119551 ) and Scalo (119981 ). 

One of the results most sensitive to the choice of the IMF is 
the metallicity gradient in the disk. In fact, the magnitude and the 
slope of chemical abundance gradients are related to the number 
of stars in each mass range, and so to the IMF. With our set of 
chemical yields (see Sect. 13.41 ), we find the best agreement with 
the observed gradients adopting a two power-law IMF analogous 
to Kroupa's IMF (with a slope -2.35 for M > 1 M Q and -1.2 
for M < 1 M ). This is the IMF adopted in this work. Other 
parameterizations, with the exception of Chabrier's IMF, result 
in general in an under-production of massive stars, and, conse- 
quently, a lower production of metals. 

We also examine the effects of a possible dependence of the 
slope of the IMF on galactocentric radius. On the basis of HST 
photometric observations of stellar clusters embedded in several 
giant Hn regions of M33, Lee et al. (120021 ) found that the IMF be- 
comes steeper (less rich in massive stars) with increasing galac- 
tocentric radius and decreasing metal abundance. Assuming a 
variable slope in the range 10 M Q < M G < 100 M as sug- 
gested by Lee et al. ( 120021 ). we find however that the resulting 



2 Recently, Hirschi et al. ( 12005b computed stellar yields for massive 
stars of solar metallicity considering also the effects of stellar rotation 
(Vrot = km s~' Vro, = 300 km s~' ), but only for a few chemical elements 
(not including S). 



metallicity gradient is not consistent with the observations, being 
characterized by a much shorter scale length than the observed 
distribution of metals. 



4. Results of the model 

In this section we present the results of our GCE model for M33 
described in Sect. [3j and compare them to the data discussed in 
Sect. [2] In particular, we examine the SFR, the distribution of 
atomic and molecular gas, and the chemical abundance gradi- 
ents. 

4.1. The evolution of the SFR 

An important constraint to GCE models is the magnitude and 
the radial profile of the SFR at the present epoch. This observ- 
able allows one to discriminate between different evolutionary 
scenarios that result in the same radial distributions of stars, gas 
and chemical abundances. 

This is illustrated by the two models presented in Fig. Q] and 
in Table [2 labeled collapse and accretion, differing in the time 
behavior of the accretion rate on the disk. In the former, the 
disk is formed by the collapse of the gas initially present in the 
galaxy's halo, and the resulting accretion rate of the disk de- 
creases rapidly with time (Fig. [TJ, simulating a rapid collapse 
phase. Since the exponential time scale of collapse is longer for 
the outer, less dense, regions of the halo, the disk is formed in 
an inside-out fashion, as clearly shown by the radial profile of 
the infall rate in Fig. [TJ>. In the accretion model, the infall rate 
is approximately constant for the entire evolution of the galaxy 
(except the inner few kpc, where it decreases by about one order 
of magnitude). In this case, the disk of M33 is built by continu- 
ous gas supply from the external medium, a process supporting 
a substantial star formation at all radii, in agreement with the 
observations discussed in Sec. l2.4l 

Both models well reproduce the present time atomic and 
molecular gas distributions, and the stellar density profile of 
M33. The differences between the two models are evident at ear- 
lier times, because the bulk of the stellar mass settles in the disk 
earlier in the collapse model than in the accretion model, due 
to the rapid gas infall and to the high star formation rate. The 
behaviour of the SFR reflects the differences in the two models 
at the present time: the actual SFR predicted by the accretion 
model is in closer agreement with the observations than the SFR 
predicted by the collapse model. In the collapse model the accre- 
tion rate and SFR vary more rapidly with time: as a consequence 
stars form earlier and the present-day SFR is at least a factor ~ 3 
lower than what can be inferred by infrared and optical observa- 
tions. Clearly, the higher infall rate on the disk at present times 
for the accretion model results in a more vigorous SFR than in 
the collapse model. But at earlier times the SFR was higher for 
the collapse model, since the larger gas reservoir accreted earlier 
by the disk was converted into stars. The present-day SFR inte- 
grated over the disk up to R — 10 kpc is 0.2 M yr~' for the col- 
lapse model and 0.5 M Q yr~' for the accretion model, but the dif- 
ference increases in the inner regions of the disk. The integral of 
the observed SFR up to R = 7 kpc, ~ 0.4 M yr" 1 (see Sect. [23]), 
is in good agreement with the accretion model results. Also the 
intermediate-age SFR derived from PNe, ~ 0.55 M yr _! from 
1 to 4 Gyr ago, and ~ 1.1 M yr~' from 1 to 8 Gyr ago, are 
in closer agreement with the accretion model which predicts 
0.7 M Q yr -1 and 0.9 M yr~' respectively. In the rest of the paper 
we adopt the accretion model to describe the evolution of M33. 
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Table 2. Parameters of the GCE models for M33 and the MW. H, H' and \i are given at the equivalent solar radius (/? Qi m33X fo 
is given at the centre. References: a, accretion (this work); b, collapse (this work); c, Molla & Diaz (2005) for a spiral galaxy of 
morphological type N = 6; d, Molla et al. ( fl996l H9971 . e Ferrini et al. (fl992l > for the MW. 
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Fig. 3. Time evolution of the stellar surface density (line types as Fig. 4. Time evolution of the surface density of diffuse gas (line 
in Fig. [2), compared with the stellar surface density determined types as in Fig. [2j, compared with the observed surface density 
by Corbelli (120031 shaded region). of atomic hydrogen (Corbeili [2003l l. 



4.2. The radial distribution of gas and stars 



We show in Figures|3]E][5J the time evolution of the surface den- 
sities of stars, atomic, and molecular gas. The diffuse gas is ac- 
creted by the disk from the external medium following an expo- 
nential law which produces a larger amount of gas in the central 
regions of the galaxy. Clouds form out of this diffuse gas, whose 
surface density decreases with time. Stars are formed at the ex- 
pense of clouds, and the cloud surface density decreases with 
time while the stellar surface density increases. The present-time 
radial distribution of diffuse and condensed gas can be compared 
with the observed radial profiles of atomic and molecular gas, re- 
spectively. The model predicts a slightly higher surface density 
of diffuse gas in the central regions and a lower cloud surface 
density than observed. This might either imply that the forma- 
tion of diffuse gas via clouds collisions is even more inefficient 
at small radii than we assumed, or that there is an enhanced for- 
mation of molecular clouds. This might be due to additional fea- 
tures in place in the central regions, such as a small bar fueling 
gas towards the center or to a small bulge with a higher metal and 
dust abundance enhancing the molecular hydrogen fraction. The 
slightly lower abundance of clouds seen at large galactic radii 
might be associated with a change of the CO-H2 conversion fac- 
tor due to, e.g. a lower excitation rate of the CO molecule. 
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Fig. 5. Time evolution of the surface density of clouds (line types 
as in Fig. [2]), compared with the observed surface density of 
molecular hydrogen (Corbelli 2003) . 



4.3. The metallicity gradients 

We now examine the evolution of the abundance gradients of 
O, Ne, S, N, and Fe in M33. In this nearby galaxy is possi- 
ble to separate the present-day metallicity gradient, outlined by 
abundances of HII regions, A-type and B-type supergiant, and 
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Cepheids, from an older metallicity gradient, outlined by PNe 
and RGB stars. The PNe trace the chemical composition of the 
ISM over a range of galactic ages between ~ 1 to ~ 8 Gyr ago, 
at least for those elements not affected by stellar evolution in the 
mass range M <8 M Q , such as O, Ne (only as a first approxima- 
tion, cf. Marigo 2001, Magrini et al. 12005 a), and S. Abundances 
in PvBG stars trace the metallicity gradient relative to an even 
distant past (>8 Gyr ago). 

In order to avoid statistical effects due to the incompleteness 
of the various samples, in the rest of this Section we shall com- 
pare the model results with larger and homogeneous samples of 
chemical abundance determinations for each class of objects, 
namely Hn regions, young stars, and RGB stars. As discussed 
by e. g Vilch ez et al. ( 119881 ), Urbaneja et al. d2005l l, Magrini 
et al. (120071) . the gradient in M33 does not have not a constant 
slope, and it cannot be represented by a single power-law: stud- 
ies undertaken for different radial ranges would produce differ- 
ent results. Inclusion of the innermost regions would result in 
a steeper gradient. On the contrary the gradient flattens out if 
only regions far from the center were considered. Thus, for each 
chemical element we compute the gradient considering all data 
from the literature within a given radial range, and we include 
also new abundance determinations by Magrini et al. (2007]). In 
Table l4.3l we list the slope of the gradients and the references for 
the data used to derive it. 



4.3.1 . The abundance of oxygen 

In Fig. [6] we compare the time evolution of the O gradient 
with observations of Hn regions (a), stars (b) and PNe (c). 
The O/H gradient predicted by the model at present time is 
-0.067 dex kpc~' from 1 to 10 kpc in radius. This is in agree- 
ment with the gradient derived in the same radial range from a 
combined sample composed by abundances in Hn regions avail- 
able in the literature and from recently acquired data presented 
in Magrini et al. ( 120071 ). A linear fit of all H n regions data, which 
references are quoted in Table l4.3l gives -0.07 ± 0.01 dex kpc -1 
over a range of radii from 1 to 10 kpc, even though the data show 
that the gradient flattens out going radially outwards. We stress 
the need to assemble a statistically significant sample in order 
to avoid the influence of intrinsic peculiarities of the observed 
sources and inhomogeneities associated to particular regions e.g. 
to spiral arms. 

The correspondence between the predicted O/H gradient and 
that determined from absorption lines in young stars is also very 
good. Monteverde et al. (|2D00l > find -0.078 ± 0.06 dex kpc" 1 , 
and a similar value was also found by Urbaneja et al. ([2005 ), 
-0.06 + 0.02 dex kpc 1 . The overall linear fit to O abundances 
from young stars, including A, B giant stars and Cepheids, (ref- 
erences quoted in Table 14.31 1. gives -0.07 + 0.02 dex kpc~' 
(Pearson correlation factor -0.6), in agreement with the model 
result (see Table |4.3l ). 

The O/H gradient of a sample of PNe reaches a lower ab- 
solute value of abundances at large radii and is steeper than 
the O/H gradient outlined by Hn regions and young stars. A 
weighted linear fit gives a slope -0.10 + 0.05 dex kpc -1 (with 
Pearson correlation factor -0.7). This behaviour is expected in 
a scenario where metallicity gradients flatten with time because 
chemical abundances at large radii increase gradually with time 
while the enrichment process is very fast in the central regions. 
The O/H gradient from PNe is representative of the metallicity 
in the disk of M33 about 1-8 Gyr ago. The same is shown by 
our model (compare the long-short dashed curve and the short 
dash-dotted line in Fig(6]:). More quantitatively, the slopes of the 
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Fig. 6. The time evolution of the O gradient (line types 
as in Fig. [2) compared with: a) Hn region determinations: 
Smith ( 119751 filled squares); Kwitter & Aller ( 119811 crosses); 
Vilchez et al. (119881 fille d circle s); Crockett et al. (120061 empty 
circles); Magrini et al. (2007] filled triangles); b) B super- 
giant stars: Monteverde et al. (1997, 2000, filled triangles); 
Urbaneja et al. (120051 empty triangles); Cepheids: Beaulieu 
et al. ( 120061 empty squares); c) PNe abundances: Magrini et 
al. d2004l crosses). 



O/H gradients predicted by our model are -0.078 dex kpc 1 and 
-0.094 dex kpc -1 , 5 and 8 Gyr ago respectively (see Table |4.3l l. 

Our sample includes some of the brightest PNe of M33, 
probably representative of a younger population than the aver- 
age sample (Richer et al. 1998 ). Five PNe show an O abundance 
in agreement with the model predictions, four are marginally 
consistent, and two are over-abundant in O (PN 75 and 91, see 
Magrini et al. 2004). A possible explanation for the O over- 
abundance in some PNe is the occurrence of hot bottom burn- 
ing and third dredge-up processes during the post-AGB phases, 
with a consequent O enrichment in the nebula (Marigo 2001 
Herwig 2004). These processes are generally associated with 
massive progenitors and therefore with N enrichment as well. 
We have therefore excluded the two PNe over-abundant in O 
from the weighted linear fit shown in Table |4.3l 

4.3.2. The abundance of neon 

O and Ne are both produced mainly by short-lived, massive stars 
(M > 10 M Q ). Stellar nucleosynthesis models (e.g. Chieffi & 
Limongi 2004 Hirschi et al. 2005), predict that the abundances 
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Table 3. Metallicity gradients (dex kpc ') from 1 to 10 kpc. The observed gradients are computed from the data discussed in 
Sect. 12.51 using a weighted linear fit. Model gradients are approximated by linear fits. References: a, Smith ( 119751 1; b, Kwitter 
& Aller ( fT98Tb : c, Vilchez et a l. ( (19881 ; d, Crock ett et al. (l20Q6t ; e, Magri ni et al. (l20QTb ; f, Will ner & Nelson-Patel (120021; g , 
Monteverde et al. ( fT997l [2000l >; h, Urbanej a et al . ( 120051 ; k, Beauli eu et a l. d20"06l ; Ma grini e t al. d20"04l ; Ba rker et al. ([2006); 
/, Stephens & Frogel d2"0"02l ; m, Kim et al. (120021 ; n, Galletti et al. (120041 ; o, Tiede et al. ( |2004] l; p, Brooks et al. (120041 . 



of these two elements should be closely correlated. This predic- 
tion is supported by abundance measurements in Galactic and 
extragalactic PNe (Henry [T990I I, showing that Ne/O is constant 
over a wide range of [O/H] values. Fig|6] and [7] show indeed 
that the slopes of the two gradients predicted by our model are 
indistinguishable and equal to -0.067 dex kpc -1 . 

Willner & Nelson-Patel (2002 ) derived the Ne abundance for 
25 Hn regions in M33 from infrared spectroscopy. They find 
an inconsistency with the O gradient measured by Vilchez et 
al. (119881 -0.12 + 0.01 dex kpc 1 . Excluding the inner two and 
the outer three regions of their sample they found a best-fit slope 
of -0.07 + 0.02 dex kpc -1 , in agreement with our predictions. 
An even shallower slope for the [Ne/H] gradient is found by 
Crockett et al. (120061 1, -0.016 dex kpc 1 . Both sets of obser- 
vations are plotted in FigfTJfor a comparison with our predicted 
gradient. Although there is some dispersion in the data in the 
central and outer regions, the overall slope of the observed Ne/H 
gradient is -0.06 + 0.02 dex kpc 1 (with Pearson correlation fac- 
tor -0.7). This has been computed using data available in refer- 
ences given in Table 14.31 and is in agreement with the gradient 
of -0.067 dex kpc 1 predicted by our model. 
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Fig. 7. The time evolution of the Ne gradient (line types as in 
Fig|2]l. Comparison with the observations: Hn regions infrared 
spectra, Willner & Nelson-Patel (2002 , filled squares) and opti- 
cal spectra, Crockett et al. ( 120061 empty circles). 



4.3.3. The abundance of sulfur 

The abundance of S in PNe is expected to be a better tracer of 
the chemical composition of the ISM at the time of formation 
of PNe progenitors than the abundance of O and Ne (Henry et 
al. l2004l l. For this reason, the S/H gradient in PNe has been used 
in our Galaxy to determine the temporal behaviour of the metal- 
licity gradient (Maciel et al. 2005 2006 ). We recall however the 
possible sulfur "anomaly" seen in some Galactic PNe (Henry et 
al . 120041 1: when compared with similar data for stars and Hn re- 
gions, some Galactic PNe have much lower S abundance than 
expected on the basis of their O abundance. This anomaly is not 
dominant in the sample of PNe by Magrini et al. (2004). The 
average log(S/0) = -1.92 is consistent with Galactic studies of 
both PNe and Hn regions (cf. -1.91 + 0.24, Henry et al. l2004l . 
It might be present at most in two PNe (PN93 and PN96) with 
very low S/O ratio. 

In FigJBJj and b we show the S abundance of Hn regions and 
PNe, respectively. Notice the good agreement of the present-day 



S gradient predicted by our model (-0.07 dex kpc -1 ), with the 
S abundances gradient for Hn regions (-0.07 + 0.05 dex kpc -1 , 
Pearson correlation factor of -0.6). A best-fit to the PN sample 
gives an S gradient of -0.09+0.04 dex kpc -1 (with a Pearson cor- 
relation factor -0.8) consistent with the value -0.082 dex kpc 1 
predicted by our model 5 Gyr ago. We have excluded the two 
PNe possibly affected by the sulfur "anomaly". These do not 
correspond to the two PNe excluded from the oxygen gradient, 
since the two phenomena, sulfur "anomaly" and oxygen over- 
abundance, are not necessarily linked. 

4.3.4. The abundance of nitrogen 

Unlike O, Ne and S, nitrogen is mainly produced by intermedi- 
ate mass stars. Therefore, the N abundance of B stars (Urbaneja 
et al. 120051 ) and PNe (Magrini et al. 120041 ) is affected by nucle- 
osynthesis processes, and cannot be used to infer the original 
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Fig. 8. The time evolution of the S gradients (line types as 
in Figf2]i compared with the observations of: a) Hn regions: 
Kwitter & Aller ( I198U cro sses); Vilchez et al. 0988, filled cir- 
cles); Magrini et al. (2007 , filled triangles); b) PNe: Magrini et 
al. (120041 crosses). 



chemical composition of the ISM. In Fig|9] we show the model 
evolution of the N/H gradient and the abundances observed in 
Hn regions. The predicted slope of the present-day N gradient, 
-0.125 dex kpc 1 , is in good agreement with Hn regions data 
(-0.10 ±0.02 dex kpc 1 , with a Pearson correlation factor -0.8). 

The general agreement of the observed gradients of chemi- 
cal elements produced by stars in different mass ranges with the 
model supports the reliability of the adopted IMF. Any choice 
of the IMF that predicts the correct number of massive stars, 
but not of low- and intermediate mass stars, would reproduce 
equally well the O/H, S/H and Ne/H gradients, but would fail in 
predicting the N radial distribution. 

4.3.5. The abundance of iron 



Barker et al. (2006) analyzed [Fe/H] abundance^] in RGB 
stars of M33 from new observations and from previous works 
by Steph ens & Frogel (2002), Kim et al. (120021 ), Galletti et 
al. ([2004] ). Tiede et al. (120041 ). Brooks et al. (120041 ). Using these 
data, they found a well defined gradient [Fe/H] = -0.07/? -0.49 
extending up to R as 12 kpc. FigJTOj shows the good agreement 
between the time evolution of Fe/H predicted by our model and 
the RGB abundances of Barker et al. (2006). The epoch of for- 
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Fig. 9. The time evolution of the N/H gradient according to 
our model (line types as in Fig. 0), compared with measure- 
ments of N in Hn regions: Smith ( 1975] filled squares); Kwitter 
& AUer ( flWTl cros ses); Vilchez et al. JT9881 filled circles); 
Magrini et al. (2007, triangles). 




3 As usual, [Fe/ H] = lo g(Fe/H)-log(Fe/H) e and 12 + log(Fe/H) e 
7.45 (Asplund et al. l2005l 
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Fig. 10. The time evolution of the [Fe/H] gradient according to 
our model (line types as in Figf2]), compared with the [Fe/H] 
gradient of RGB stars (Barker et al. 120061 shaded area). 



mation of RGB stars is >8 Gyr ago, when the slope of the Fe/H 
gradient according to our model is -0.08 dex kpc -1 . Note how- 
ever that the linear fit is a very rough approximation to the effec- 
tive model gradient. 

The metallicity gradient of Cepheids determined by Beaulieu 
et al. 120061 -0.2 dex kpc" 1 , can also be compared with the 
predicted gradient at present time. However, the five observed 
Cepheids are located in the radial region from ~1 to 4 kpc, 
whereas our computed gradient extends out to 10 kpc. In the 
inner region of M33, the [Fe/H] gradient predicted by our ac- 
cretion and collapse models model are -0.11 dex kpc" 1 and 
-0.14 dex kpc" 1 , respectively. We must keep in mind however 
that the gradient quoted by Beaulieu et al. (2006) refers to the 
total metallicity, not to the Fe/H gradient. 
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4.4. Comparison with previous models 

Specific models for the chemical evolution of M33 have been 
previo usly developed by Dia z & To si d!984j ), Molla et al. d 19961 
119971 , and by Molla & Diaz ([2005 ). Our approach is very similar 
to that of Molla et al. d 19961 [1997 ). and the predicted H I and H 2 
radial distributions are indeed similar. The chemical gradients 
predicted by the model of Molla et al. ( 119961 1 19971 1 are how- 
ever steeper (e.g. d(0/H)/dR ~ -0.21 dex/kpc, d(N/H)/dR ~ -0.34 
dex/kpc) than what the most recent abundance determinations 
suggest. This is probably a consequence of the assumed IMF, 
heavily weighted towards massive stars in the parametrization 
adopted by Molla et al. d 1 9961 II 9971 . The time evolutio n of th e 
metallicity gradients discussed in detail in Molla et al. d 19971 1, 
shows the same behaviour as in the present work, i.e. a flat- 
tening with time. Metallicity gradients resulting from the most 
recent models by Molla & Diaz (2005), not computed specif- 
ically for M 33 but applied to the M 33 case, are also steeper 
than what the collection of observations presented here seem to 
show. Molla & Diaz (|20051 > assume a slightly higher infall rate 
than Molla et al. ([1996 1997 ), resulting in significantly higher 
surface density of H I and H2, which are marginally consistent 
with the data. In particular their model predicts a strong decrease 
of the H I surface density towards the center, a large central H I 
hole, and a slight decrease of the H 2 surface density in the same 
region, which are not observed in this galaxy. 



5. The time evolution of metallicity gradients 

Both models we consider predict abundance gradients which 
flatten with time. The main difference between the two model 
is the resulting value of the slope of the gradients, that are much 
steeper at all times in the collapse model than in the accretion 
model. This is due to the different nature of the infall: the col- 
lapse model predicts a much shorter time scale for the formation 
of the disk, due to a rapid collapse of the halo, while the ac- 
cretion model predicts a continuous infall of material from the 
intergalactic medium. 

In the case of M33, the best fit to the whole set of observa- 
tions is obtained with our accretion model, or, in the terminology 
by Molla et al. 119931 with a model which has a self consistent 
SFR and a continuous gas infall. This model does not exclude 
an early collapse phase of a galactic halo, but considers its ef- 
fects less important at the present time than the accumulation of 
matter from the intergalactic medium. This is supported by re- 
cent observations (see Sect. |2.4| i and by numerical simulation of 
galaxy formation. As discusses in Sect. 13. li the intergalactic sup- 
ply is important, especially in low mass halos, and, in the case 
of M33. The continuous supply of gas is necessary to reproduce 
the high SFR observed at present time (see Fig. [2j. 

In Fig. [TT] we compare the time evolution of the slope of 
the gradient of O/H and S/H (panel a, accretion model; panel b, 
collapse model) with the gradient measured in young stars, Hn 
regions and PNe. To derive the gradient, both the model predic- 
tions and the observational data have been linearly fitted in the 
log[M I H]-\og R plane, in the range of galactic radii 1-10 kpc. 
The flattening of the metallicity gradient predicted by the accre- 
tion model over the entire lifetime of the disk of M33 is quite 
modest because of the slow process of formation of the disk, 
which is still taking place at the present time. On the other hand, 
the collapse model predicts steeper gradients at all times, espe- 
cially in the past, because of the rapid initial collapse phase that 
forms the disk in a few Gyr. This behaviour is in conflict with 
the observations, especially with the Fe/H gradient of RGB stars 
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Fig. 11. The time evolution of the [O/H] (a), [S/H] (b), and 
[Fe/H](c) gradients according to our accretion (continuous line) 
and collapse (dotted line) models, compared with gradient 
slopes of sulfur (filled triangle) and oxygen (filled square) of 
PNe, of sulfur (empty triangle) and oxygen (empty square) of 
Hn regions and of young giant stars (star), and of iron of RGB 
stars (cross). The vertical dashed line marks the present age of 
the galaxy, 13.6 Gyr. 



(see Sect. |4.3.5l l determined by Barker et al (2006), as shown in 
Fig. [CD:. 

The comparison of the abundances in PNe with that in Hn 
regions presented in Sect. 14.31 represents the first evidence for 
a flattening of the abundance gradients with time in an external 
galaxy. A similar result has been found in the MW by Maciel 
et al. d2006l ). The flattening in M33 is particularly evident for 
O and S abundances, as shown by Fig. [6] and Fig. [8] Another 
feature of particular interest in the case of M33 is the overall 
increase of metallicity with time at all radii. The steady infall 
rate, that increases the gas content of the disk with time, fuels 
star formation at all radii. The flattening of the gradient with time 
is due to the combination of the weak radial dependence of the 
infall rate and to the star formation efficiency which decreases 
more rapidly with time in the central regions. 



6. Conclusions 

We have analyzed a large sample of observations for the LG 
spiral galaxy M33, including gas and stellar radial distributions, 
metallicity gradients, and star formation rate, and we have inter- 
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preted them using GCE models for this galaxy. We summarize 
here our main results: 

(z) Our accretion model (where the infall rate on the disk is al- 
most constant in time) reproduces the data better than our 
collapse model (where the accretion rate decreases rapidly 
to simulate the early collapse phase of a baryonic halo). 
Whereas both models can account (within the errors) for the 
observed distribution of stars, atomic and molecular gas, the 
abundance gradients, outlined by a large number of abun- 
dance determinations, and the high value of the present-day 
SFR allow one to discriminate between different accretion 
and SFR histories for this galaxy. 

(z'z) The ability of our accretion model to reproduce the observed 
constraints suggests the existence of an extended phase for 
the formation of the disk of M33. The continuous accretion 
of material from the intergalactic medium, is also supported 
by deep observations at 21-cm in the proximity of M33 and 
by numerical simulation of galaxy formation and evolution. 
The continuous infall of gas into the disk is necessary to 
maintain a high SFR, declining very slowly with time. 

(Hi) Chemical abundances determined in young stars, Hn regions, 
PNe, and RBG stars, all indicate a relatively fiat radial distri- 
bution of elements, suggesting a slow, continuous rate of star 
formation in the disk of M33. The magnitude and evolution 
of the gradients of O/H, Ne/H, N/H, S/H, [Fe/H] computed 
with our accretion model are in good agreement with these 
observations. In particular, S/H and O/H gradients of PNe 
when compared with the same elemental abundance gradi- 
ents of HII, and the [Fe/H] gradient of RGB stars, give inde- 
pendent observational evidence that the metallicity gradient 
has flattened over the last ~ 8 Gyr, in agreement with our 
model results. A similar behaviour has been recently found 
in the MW (Maciel et al. 2006) and may represent a general 
feature of the evolution of disk galaxies, even though galax- 
ies of different masses, such as the MW and M33, might have 
followed different evolutionary histories. 
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Appendix A: Chemical abundances in M33 

We present here the collection of chemical abundance measure- 
ments used in the present work. 
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Table 1. Chemical abundances used to constrain the present chemical evolution model of M 33 from H II regions (S75, Smith [19751 KW81, 
Kwitter & Aller HMTl V8 8, Vilchez et al. [T988l C06 , Crockett et al. [20061 M0 7, Mag rini et al . 120071 . young stars (M97, Monteverde et al. [T997l 
U05, Urbaneja et al. 2005 ; B06, Beaulieu et al. 2006 1, PNe (M04, Magrini et al. 2004 1. Chemical abundance marked with * are from recompilation 
by Garnett et al. ( 1997) and those marked with : are derived without the electron temperature measurement. We adopted typical errors of 0.2 dex 
when they are not quoted. 
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+jU 4 / Uj.o 


8.35* 


7.30* 






S75 




Mpp coo 


1 JZ 4j.Z 


_i_"iPi ISl C/l Pi 
+jU jo J4.U 


o.4j 


7.00* 






S75 




IC 132 


1 33 15.8 


i 1C\ CA /I C A 


8.05* 


6.75* 






S75 




dot f \ x K ii\ 
BLLMr oj (MA 11) 


T 1 A f\n A 

1 34 0/.0 


i "3 A H i T3 

+30 3y Z3 


8.70* 


7.35* 




7.1 1 


KW81 : 




IC 142 


1 33 55 1 




8 72* 


7 50* 


7 42 


6 78 


KW81 ■ 




INLiC jVj 


1 33 35.50 


+30 41 52 


8.22 : - : 


7.00* 


7.75 


7.0 1 


KW81 : 




dLLMF «o (MA 2) 


1 34 15.5 


+30 37 1 1 


8.75* 


7.50* 




6.78 


KW81 : 




MA 3 


1 34 01 


+30 52. 1 


8.45 


7.10 


7.63 


6.75 


KW8 1 




Mpp AA/I 
INOC oU4 


i i/i 1 1 in 

1 J4 jj.iy 


i 1C\ A "7 A< A 


8.30 


6.97 


7.49 


6.70 


KW8 1 




IC 131 


1 33 15.0 


i an /i c ao 
+JU 4j uy 


8.05 


6.72 


7.09 


6.60 


KW8 1 




Dpi \/IT> ACA /M A Q-i\ 

DLLrvir oju (iviA yaj 


1 1A 11 f^A 
i J4 3 J. 04 


i Q i AA 1 T 
+3 I UU Zl .z 


Q CA* 

o.jU 




7 11 


7.10 


KW8 1 : 




IC 133 


1 33 15.9 


i "JA C'Q AT 
+JU jj uz 


H.25'- : 


6.9* 


7.55 


6.95 


KW81 : 




Mpp coo 
INLiC joo 


1 32 45.2 


i 1l\ 1Q C/1 A 
+JU jo J4.U 


S.( )(.) 


6.84 


7.30 


6.68 


KW8 1 




tp i n 

1L. 1 JZ 


1 11 K Q 
1 JJ iJ.o 


, Ifl Cfl AC a 

+JU JO 4j.U 


Q AA 


O. 1 z 


7 1Q 


A 7A 
O. I\) 


KW81 




BLLMr 93 (CCy3) 


1 33 52.3 


■ n\ in in 
+30 39 lo 


y.oz±o.i6 


7.92 




n 3 1 i A Hi" 

/. 23 ±0.06 


V88: 




IC 142 


1 33 55.1 


+30 45 22 


8.70±0.16 


7.57 




n Al i A AC 

7.03±0.05 


V88 




NGC 595 


1 1 -3 -3 C CA 

1 33 35.50 


+30 41 52 


O A A i A AA 

8. 44 ±0.09 


7.33 




£. OZT i A AO 

6.86±0.0a 


V88 




BCLMP 88 (MA 2) 


1 34 15.5 


+30 37 11 


8.44±0.15 


7.24 




6.8±0. 1 


V88 




NGC 604 


i iii n in 

1 34 33.19 


i -3 A A 1 AC H 

+30 47 05.6 


O C 1 i A A "3 

8.5 1±0.03 


7.35 




£ AC i A A-3 

6.95±0.03 


V88 




IC 131 


1 33 15.0 


+30 45 09 


8.41±0.06 


7.25 






V88 




NGC 588 


1 32 45.2 


+30 38 54.0 


8.30±0.06 


6.77 




6.86±0.()6 


V88 




BCLMP 027 


1 33 45.5 


30 36 5 1 






7.68±0.09 




W03 




BCLMP 079 


1 34 0.2 


30 40 49 






7.68±0.04 




W03 




BCLMP 087E 


1 34 2.5 


30 38 41 






7.96±0.08 




W03 




BCLMP 042 


1 33 5.5 


30 39 30 




- 


7.89±0.10 


- 


W03 




BCLMP 301 


1 33 5.3 


30 45 22 






7.63±0.08 




W03 




BCLMP 004 


1 33 9.3 


30 35 49 


■ 


- 


7.85±0.08 




W03 




BCLMP 062 


1 33 4.5 


30 44 38 


■ 


- 


7.55±0.08 




W03 




BCLMP 049 


1 33 3.9 


30 41 28 


■ 


■ 


7.77±0.04 




W03 




BCLMP 045 


1 33 8.8 


30 40 25 


"" 




7.66±0.04 




W03 




Dpi \/lT> 
BLLlVlr JUZ 


i o4 o.y 


n\ An 0£ 

3\) 4 / Zu 






Q AA-lA A1 




W03 




BCLMP 214 


1 33 30 


30 3 1 47 






-7 O A i A no 

7.84±0.08 




W03 




BCLMP 095 


1 34 0.9 


30 36 18 






n on i a no 
/.Q/+U.UQ 




W03 




BCLMP 088W 


1 34 5.5 


30 37 12 






7.72±0.05 




W03 




BCLMP 710 


1 34 3.6 


-3 A -3 -3 A ^ 

30 33 42 






7.87±0.07 




W03 




BCLMP l\Jl 


1 34 10 


30 31 57 






h oa ■ n no 
/.o4±0.0o 




W03 




RPI MP AO, 1 


1 34 6 4 


30 51 55 






7 5K+0. OA 




W03 




Dpi \/IT> ACAP 
DL.L,lVlr OoUL 


1 J4 Z. 1 


1f\ A "7 A 
3\) 4 / U 






7 7A-1-A (M 




W03 




DPT NflD A^IID 

BLLMr OoUd 


1 34 3.5 


30 46 50 






7 AO_i_A Al 

/.oZ±U.Uj 




W03 




DPI MD T) 1 

BLLMr ZZ1 


1 33 9.8 


in -1-7 <-ic 

jU Zi Zj 






7 AC_i_A A7 
/.OJ±U.U/ 




W03 




DPT HiTD7/1AM17 

BLLMr /4UW 


1 34 9.6 


30 41 52 






7 A A_i_A AC 
/.4o±U.Uj 




W03 




LrJJor iy4 


1 33 1.1 


1 A /I c, 1 £ 
JU 4j 1 






7 Ci-i-A A7 
/.JJ±U.U/ 




WO 3 




DPI MD TCI 

BLLMr Zj 1 


1 3336.6 


in oa i i 
3U ZU 13 






A AA_i_A 1 A 
O.UU±U. 1U 




W03 




DPT KlTD^ll 

BLLMr OZJ 


1 3316.5 


in c3 cn 

jU jZ j\J 






A AA_i_A Al 

o.yu±u.uj 




W03 




dpi ivad ogA 
BLLMr ZoU 


1 3245.2 


in ig c/i 

jo j4 






7 1A_i_A AA 




WO 3 




BCLMP 638b 


1 33 16.3 


30 56 44 






n a i i n no 

/.41±0.0s 




W03 




Dpi \/lD ATflM 


1 33 1 J.D 


1A CA /1Q 

ju jo 4y 






7 CA-lA AS 
/.jU±U.US 




W0 3 




BCLMP 090 


1 34 04.2 


i "> n 1 O A A -3 

+30 38 09.2 


8.50±0.06 




7.96±0.09 




C06 




BCLMP o91 


1 34 16.6 


i 1 A C 1 C A t \ 

+30 51 54.0 


O -tiC i A A-t 

o.26±0.02 




1 tin i a m 
/.62±0.03 




C06 




BCLMP 745 


1 34 37.6 


i -3 A 1 /I CC A 

+30 34 55.0 


8.07±0.10 




7.2±0.2 




C06 




RPT MP 7AA 

DV^LilVir / \J\J 


1 34 42 2 


j. ■an 31 d.9 % 

i 3\) J 1 H-Z, . J 


Q.JZiU. 1 Z 




7 6±0 3 




C06 




Dpi \/\r> OQA 
tSL-LiVlr ZyU 


1 11 1 1 A 

l 33 1 1 .4 


,1(1 /If IC 1 

+JU 4j 1 J. 1 


Q 1 J.A AA 
o.Zl±U.UD 




7.4±0. 1 




PAA 
L UO 




MAI 


i 11 Al A 
I 33 UJ.4 


+jU 11 10./ 


O -I A ,A Ait 

o.Z4±U.Uo 




7 A I _i_A AQ 




C06 




LGC HII 2 


1 -3 -3 A -3 n 

1 32 43.0 


1 A 1 A -3 1 -3 

30 19 31.2 


O -3 C i A (\H 

8.25 ±0.06 


7.17±0.25 




6.75±0. 18 


M07 




nil J 


1 32 45 9 


30 41 35 5 


O.itlU.UJ 


A Q4-t-A 1 7 




A A9-1-A 1 1 

O.OZiU. 1 1 


M07 




Dpi 


1 IT CO C 

1 JZ JfiJ 


1A AA OQ (. 
JU 44 Zo.O 


O.Zj±U. Yd 


A OAj-A c,a 

o.yu±u. JU 




A Qc,A 
O.o J±U. J 


M07 




DPT Hif DO 1 O 

BCLMrZl 


1 "31 Aft T 


ic\ n\ An i 
jU jU 4 / .j 


g -3 c_i_a AC 
S.Zj±U.Uj 


n i n i a i a 

i . iy±u. 14 




6.77±0. 15 


M07 




ivi l, ivi uu t mz4 




1A 1 C AQ c, 
JU 1 o Uo.J 


C 1 Cj-A Oc, 
o. lo±U.ZJ 


7 1 A j_A 1A 
/. 1U±U.jU 




A CC-lA 1 
O.JjiU. J 


M07 




CPSDP194 


1 33 11.1 


30 27 34.2 


8. 27 ±0.06 


7.12±0.20 




6.67±0. 17 


M07 




BCLMP626 


1 33 16.4 


30 54 04.8 


8.13±0.07 


6.86±0.30 




6.61±0.28 


M07 




BCLMP45 


1 33 29.0 


30 40 24.8 


8.49±0.04 


7.17±0.15 




6.92±0.17 


M07 




BCLMP637 


1 33 50.6 


30 56 33.3 


8.34±0.05 


7.30±0.40 




6.90±0.40 


M07 




GDK99 128 


1 33 59.9 


30 32 44.3 


8.47±0.06 


7.31±0.21 




6.98±0.21 


M07 




BCLMP670 


1 34 03.3 


30 53 09.3 


8.28±0.07 


7.14±0.30 




6.75±0.18 


M07 




VGHC 2-84 


1 34 06.7 


30 48 56.4 


8.35±0.04 


6.95±0.18 




6.56±0.12 


M07 




BCLMP717b 


1 34 37.4 


30 34 54.3 


8.18±0.07 


7.04±0.28 




6.57±0.24 


M07 




LGC HII 1 1 


1 34 42.2 


30 24 00.5 


8.17±0.06 


6.68±0.26 




6.42±0.40 


M07 


giant stars 


M33 1054 


1 33 50.61 


30 38 36.70 


9.03 








M97 




M33 1345 


1 33 59.676 


30 23 00.44 


8.50 








M97 




M33 B133 


1 33 28.848 


30 47 46.32 


8.56 








M97 




M33 110A 


1 33 41.023 


30 22 36.97 


8.43 








M97 




10900 


1 33 44.90 


+30 36 16.70 


8.30±0.40 


8.75±0.36 






U05 




1110A 


1 33 41.00 


+30 22 37.00 


8.40±0.40 


8.50±0.42 






U05 




1B38 


1 33 00.83 


+30 35 05.10 


8.15±0.29 


8.50±0.20 






U05 




1B133 


1 33 29.00 


+30 47 44.00 


8.10±0.16 


8.3()±0.17 






U05 
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Table 1. continued. 



Magrini et al.: The evolution of M33 



Type 


Name 


RA 


Dec 


O/H 


N/H Ne/H 


S/H 


Ref. 






J2000.0 












1 1054 


1 33 50.83 


+30 38 34.50 


8.10±0.15 


8.90+0.18 - 


- 


U05 




11137 


1 33 53.23 


+30 35 26.10 


8.50±0.30 


8.70±0.25 - 


- 


U05 




10B11241 


1 33 42.02 


+30 21 42.30 


8.00±0.20 


8.60+0.15 


- 


U05 




1 UIT1 03 


1 33 27.29 


+31 00 56.70 


7.90±0.10 


8.20+0.13 


- 


U05 




l urn 22 


1 33 33.70 


+30 47 20.20 


8.10±0.12 


8.50+0.06 


- 


U05 




1UIT136 


1 33 35.73 


+31 00 47.00 


7.70±0.19 


8.40+0.18 


- 


U05 




1OB10-10 


1 33 44.22 


+30 31 48.20 


8.55+0.21 


8.65+0.22 - 


- 


U05 


Cepheids 


A121.029 


1 34 59.72 


30 52 25 2 


8.22 


- 


- 


B06 




B160.520 


1 32 56 82 


30 41 33 8 


8.30 


- 


- 


B06 




C133.449 


1 34 33 43 


30 51 1 5 fi 


8.56 


- 


- 


B06 




D234.922 


1 33 S4 63 


30 35 1 Q 8 


8.57 


- 


- 


B06 




E237.367 


1 34 03.97 


30 38 08.4 


8.58 


- 


- 


B06 


PNe 


PN8 


1 32 28.66 


30 25 53.2 


7.94 


7.07 


6.38 


M04 




PN18 


1 33 06.11 


30 31 04.5 


8.20 


7.48 


6.83 


M04 




PN28 


1 33 19.25 


30 29 40.4 


8.69 


6.98 


6.47 


M04 




PN60 


1 33 46.42 


30 26 55.3 


8.69 


7.67 


6.57 


M04 




PN65 


1 33 49.37 


30 32 06.3 


8.67 


7.61 


6.84 


M04 




PN75 


1 34 01.14 


30 50 27.1 


8.49 


8.03 


6.81 


M04 




PN91 


1 34 13.97 


30 22 36.4 


8.48 


8.07 


6.60 


M04 




PN93 


1 34 15.88 


30 24 54.6 


8.47 


7.75 


6.66 


M04 




PN96 


1 34 15.44 


30 32 20.2 


8.76 


7.43 


6.44 


M04 




PN101 


1 34 24.25 


30 27 54.3 


8.72 


7.33 


6.18 


M04 




PN125 


1 34 32.85 


30 41 10.5 


8.31 


7.53 


6.66 


M04 



